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Abstract 

We propose a model of fractal point process driven by the nonlinear stochastic differential equa- 
tion. The model is adjusted to the empirical data of trading activity in financial markets. This 
reproduces the probability distribution function and power spectral density of trading activity ob- 
served in the stock markets. We present a simple stochastic relation between the trading activity 
and return, which enables us to reproduce long-range memory statistical properties of volatility by 
numerical calculations based on the proposed fractal point process. 
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I. INTRODUCTION 



There are empirical evidences that the trading activity, the trading volume and the 
volatility of the financial markets are stochastic variables with the power-law probabil- 
ity distribution function (PDF) □ D and the lone-range correlations SHU Most of 
proposed models apply generic multiplicative noise responsible for the power-law probabil- 
ity distribution function (PDF), whereas the long-range memory aspect is not accounted 
in the widespread models (|.The additive-multiplicative stochastic models of the financial 
mean-reverting processes provide rich spectrum of shapes for PDF, depending on the model 
parameters []|, however, do not describe the long-memory features. Empirical analysis con- 
firms that the long-range correlations in volatility arise due to those of trading activity j^J. 
On the other hand, trading activity is a financial variable dependant on the one stochastic 
process, i.e.interevent time between successive market trades. Therefore, it can be modeled 
as event flow of the stochastic point process. 

Recently, we investigated analytically and numerically the properties of the stochastic 
implicative point processes, derived formula for the power spectrum flQ aud related the 
model with the multiplicative stochastic differential equations [lflj. Preliminary comparison 
of the model with the empirical data of the power spectrum and probability distribution of 



stock market trading activity ll| stimulated us to work on the more detailed definition of 
the model. Here we present the stochastic model of the trading activity with the long-range 
correlations and investigate its connection to the stochastic modeling of the volatility. The 
proposed stochastic nonlinear differential equations reproduce the power spectrum and PDF 
of the trading activity in the financial markets, describe the stochastic interevent time as 
the fractal-based point process and can be applicable for modeling of the volatility with the 
long-range autocorrelation. 



II. MODELING FRACTAL POINT PROCESS BY THE NONLINEAR STOCHAS- 
TIC DIFFERENTIAL EQUATION 

Trades in financial markets occur at discrete times t\, t2, t\, ... and can be considered 
as identical point events. Such point process is stochastic and totaly defined by the stochastic 
interevent time = tk+i — A fractal stochastic point process results, when at least two 
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statistics exhibit the power-law scaling, indicating that represented phenomena contains 



clusters of events over all scales of time 



12j . The dimension of the fractal point process is 



a measure of the clustering of the events within the process and by the definition coincides 
with the exponent of the power spectral density of the flow of events. 

We can model the trading activity in financial markets by the fractal point process as 
its empirical PDF and the power spectral density exhibit the power-law scaling 0, In 
this paper we consider the flow of trades in financial markets as Poisson process driven by 
the multiplicative stochastic equation. First of all we define the stochastic rate n = 1/r of 
event flow by continuous stochastic differential equation 



dr = 7 r 2/ ^ 2 dt + CTT^ 1/2 diy, 



(1) 



where W is a standard random Wiener process, a denotes the standard deviation of the 
white noise, 7 < 1 is a coefficient of the nonlinear damping and \x defines the power of 
noise multiplicativity. The diffusion of r should be restricted at least from the side of high 
values. Therefore we introduce an additional term — ycr 2 f r 2/1_2 into the Eq. (JTJ), 
which produces the exponential diffusion reversion in equation 



dr 



m 



(2) 



where m and tq are the power and value of the diffusion reversion, respectively. The asso- 
ciated Fokker-Plank equation with the zero flow gives the simple stationary PDF 



P(t 



a+l 

~ r exp 



(3) 



with a = 2( 7(J — a) and 7cr = 7/cr 2 . Eq. (J2J describes continuous stochastic variable r, defines 
the rate n = 1/r and, after the Ito transform of variable, results in stochastic differential 
equation 

dn = a 2 [(1 - 7 J + - f^V] n^dt + anMW, (4) 
2 V n J . 

where 77 = ~ — // and no = 1/tq. Eq. describes stochastic process n with PDF 

P(n)~lex P {-(^) m }, \ = 2( V -l+ la ), (5) 
and power spectrum S(f) P. Q lic| 

3-2 7ct 



(6) 



Noteworthy, that in the proposed model only two parameters, 7 CT and rj (or /i), define 
exponents A and (3 of two power-law statistics, i.e., of PDF and of the power spectrum. Time 
scaling parameter a 2 in Eq. can be omitted adjusting the time scale. Here we define 
the fractal point process driven by the stochastic differential equation (@J) or equivalently 
by Eq. (J2J), i.e., we assume r(t) as slowly diffusing mean interevent time of Poisson process 
with the stochastic rate n. This should produce the fractal point process with the statistical 
properties defined by Eqs. (0) and (JHJ). Within this assumption the conditional probability 
of interevent time r p in the modulated Poisson point process with the stochastic rate 1/r is 

1 



<p{t p \t) 



cxp 



r J 



(7) 



Then the long time distribution (p(r p ) of interevent time r p has the integral form 



V(Tp) = C 



poo 


_v 




r / \ m-i 


/ exp 




T a exp 
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/o 


T - 







dr, 



(8) 



with C defined from the normalization, J °° <f(r p )dT p = 1. In the case of pure exponential 
diffusion reversion, m — 1, PDF (jHj) has a simple form 



<p(t p ) 



l + q 

'p 



K (l+a) 2 



(9) 



r(2 + a)r \t 0j 

where K a (z) denotes the modified Bessel function of the second kind. For m > 1 more 
complicated structures of distribution y(r p ) expressed in terms of hypergeometric functions 
arise. 



III. ADJUSTMENT OF THE MODEL TO THE EMPIRICAL DATA 

We will investigate how the proposed modulated Poisson stochastic point process can be 
adjusted to the empirical trading activity, defined as number of transactions in the selected 
time window r<j. Stochastic variable n denotes the number of events per unit time interval. 
One has to integrate the stochastic signal Eq. (J3J) in the time interval to get the number 
of events in the selected time window. In this paper we denote the integrated number of 
events as 

n(t')dt' (10) 
and call it the trading activity in the case of the financial market. 
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Detrended fluctuation analysis |13( is one of the methods to analyze the second order 
statistics related to the autocorrelation of trading activity. The exponents v of the detrended 
fluctuation analysis obtained by fits for each of the 1000 US stocks show a relatively narrow 
spread of v around the mean value v = 0.85 ± 0.01 Q]. We use relation /3 — 2v — 1 between 
he exponents v of detrended fluctuation analysis and the exponents (3 of the power spectrum 
"l\ and in this way define the empirical value of the exponent for the power spectral density 
P = 0.7. 

Our analysis of the Lithuanian stock exchange data confirmed that the power spectrum of 
trading activity is the same for various liquid stocks even for the emerging markets The 
histogram of exponents obtained by fits to the cumulative distributions of trading activites 
of 1000 US stocks [3] gives the value of exponent A = 4.4 ± 0.05 describing the power-law 
behavior of the trading activity. Empirical values of /3 — 0.7 and A = 4.4 confirm that the 
time series of the trading activity in real markets are fractal with the power law statistics. 
Time series generated by stochastic process are fractal in the same sense. 

Nevertheless, we face serious complications trying to adjust model parameters to the 
empirical data of the financial markets. For the pure multiplicative model, when /x = 1 or 
i] = 3/2, we have to take 7cr = 0.85 to get (3 = 0.7 and 7cr = 1.7 to get A = 4.4, i.e. it is 
impossible to reproduce the empirical PDF and power spectrum with the same relaxation 
pa— 7 . and exponent of —at ivity We have proposed possible solution of this 



problem in our previous publications [2 



111 ] deriving PDF for the trading activity iV 



P(N) ~ { " a+a ' ' ' (11) 

iV» T - 1 . 

When N ^> 7 _1 this yields exactly the required value of A = 5 + 2a = 4.4 and (3 = 0.7 
for 7(T = 0.85. 

Nevertheless, we cannot accept this as the sufficiently accurate model of the trading 
activity since the empirical power law distribution is achieved only for very high values of 
the trading activity. Probably this reveals the mechanism how the power law distribution 
converges to normal distribution through the growing values of the exponent, but empirically 
observed power law distribution in wide area of N values cannot be reproduced. Let us notice 
here that the desirable power law distribution of the trading activity with the exponent 
A = 4.4 may be generated by the model with r] = 5/2 and 7(T = 0.7. Moreover, only the 
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smallest values of r or high values of n contribute to the power spectral density of trading 

n 

activity [10]. This suggests us to combine the stochastic process with two values of \i: (i) 
fi ~ for the main area of r and n diffusion and (ii) \x = 1 for the lowest values of r or 
highest values of n. Therefore, we introduce a new stochastic differential equation for n 
combining two powers of the multiplicative noise, 



dn = a 



m /no 
7ct) + 2 In 



n 4 . cm 5 / 2 



, x9 dt + - -dW, (12) 

(ne + 1) 2 (ne+1) v ' 

where a new parameter e defines crossover between two areas of n diffusion. The corre- 
sponding iterative equation for in such a case is 



Tfc+l = T fc + 



m 2 / r 



7 cr , 

2 Vr 



+ o- Tk gfc, (13) 



(e + r fc ) 2 e + r fe 

where denotes uncorrelated normally distributed random variable with the zero expecta- 
tion and unit variance. 

Eqs. flT2j) and (JT3*Jl define related stochastic variables n and r, respectively, and they 
should reproduce the long-range statistical properties of the trading activity and of waiting 
time in the financial markets. We verify this by the numerical calculations. In figure Q we 
oresent the power spectral density calculated for the equivalent processes (|T2j) and (fEfy (see 
l| for details of calculations). This approach reveals the structure of the power spectral 
density in wide range of frequencies and shows that the model exhibits not one but rather 
two separate power laws with the exponents f3\ = 0.33 and $2 = 0.72. From many numerical 
calculations performed with the multiplicative point processes we can conclude that com- 
bination of two power laws of spectral density arise only when the multiplicative noise is a 
crossover of two power laws as in Eqs. ()12|) and (|13j). We will show in the next section that 
this may serve as an explanation of two exponents of the power spectrum in the empirical 
data of volatility for S&P 500 companies |15|. 

Empirical data of the trading activity statistics should be modeled by the integrated flow 
of events iV defined in the time interval t& ^> tq. In figure El we demonstrate the probability 
distribution functions P(N) and its cumulative form P > (n) calculated from the histogram 
of N generated by Eq. ()13|) with the selected time interval Td = 100. This illustrates that 
the model distribution of the integrated signal N Jias the power-law form with the same 
exponent A = 4.4 as observed in empirical data 
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FIG. 1: Power spectral density S(f) averaged over 100 realisations of series with 1000000 iterations 
and parameters 7 = 0.0004; a = 0.025; e = 0.07; tq = 1; m = 6. Straight lines approximate power 
spectrum S ~ l/f^- 2 with ft = 0.33 and /3 2 = 0.72: a) £*(/) of the flow I(t) = *£ k ~ *fe) with 
the interevent time rfc = ifc+i — ifc generated by Eq. Q13|). b) calculated by the Fast Fourier 
Transform of n series generated by Eq. (|12|). 

P P cumulative 
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FIG. 2: a) PDF P(N) calculated from the histogram of N generated by Eq. (J13|) with the selected 
time interval = 100. b) cumulative distribution P > (N). Other parameters are as in figure ^ 

The power spectrum of the trading activity iV has the same exponent (3 = 0.7 as power 
spectrum of n in the low frequency area for all values of r<j. 

The same numerical results can be reproduced by continuous stochastic differential equa- 
tion (fT2)) or iteration equation (JT3J). One can consider the discrete iterative equation for the 
interevent time (JT3j) as a method to solve numerically continuous equation 



dr 



m-i 



7 o" , 

2 Vt 



-J— -2 dt + a-^L-dW. (14) 



(e + r) z e + r 

The continuous equation (|12jl follows from the Eq. (|14jl after change of variables n = 1/r. 
We can conclude that the long-range memory properties of the trading activity in the 




FIG. 3: PDF of interevent time P(t): open circles, calculated from the histogram of r p generated 
by Eq. ((7J) with the mean interevent time calculated from Eq. ()14|) ; open squares, calculated form 
the iterative equation (|13f) . Used parameters are as in figure 1. Straight line approximates power 
law P(r p ) ~ r-°- 2 . 

financial markets as well as the PDF can be modeled by the continuous stochastic differential 
equation (|12p. In this model the exponents of the power spectral density, (3, and of PDF, 
A, are denned by one parameter 7^ = 7/cr 2 . We consider the continuous equation of the 
mean interevent time r as a model of slowly varying stochastic rate 1/r in the modulated 
Poisson process (JJJ). In figure El we demonstrate the probability distribution functions P(r p ) 
calculated from the histogram of r p generated by Eq. (J7J) with the diffusing mean interevent 
time calculated from Eq. ()14j) . 

Numerical results show good qualitative agreement with the empirical data of interevent 
time probability distribution measured from few years series of U.S. stock data [3]. This 
enables us to conclude that the proposed stochastic model captures the main statistical 
properties including PDF and the long-range correlation of the trading activity in the fi- 
nancial markets. Furthermore, in the next section we will show that this may serve as 
a background statistical model responsible for the statistics of return volatility in widely 
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accepted geometric Brownian motion (GBM) of the financial asset prices. 



IV. MODELING LONG-RANGE MEMORY VOLATILITY 

The basic quantities studied for the individual stocks are price p(t) and return 

x(t, r d ) = \np(t + r d ) - \np(t) (15) 

Let us express return x(t, r d ) over a time interval r d through the subsequent changes Sxi 
due to the trades i = 1, 2....N(t, r d ) in the time interval [t, t + r d ], 

N(t,T d ) 

x(t,r d ) = ^2 Sx i- ( 16 ) 
i=i 

We denote the variance of 5xi calculated over the time interval r d as W 2 (t, r d ). If 5xi are 
mutually independent one can apply the central limit theorem to sum This implies 

that for the fixed variance W 2 (t, r d ) return x(t, r d ) is a normally distributed random variable 
with the variance W 2 (t, T d )N(t, r d ) 

x(t,T d ) = W(t,T d )y/N(t,T d )£ t , (17) 

where St is the normally distributed random variable with the zero expectation and unit 
variance. 

n 

Empirical test of conditional probability P(x(t, T d )\W(t, r d )) |4j confirms its Gaussian 
form, and the unconditional distribution P(x(t, r d )) is a power-law with the cumulative 
exponent 3. This implies that the power-law tails of returns are largely due to those of 
V (t, ta ) . Here we refer to the theory of price diffusion as a mechanistic random process 
12lllfl|- For this idealized model the short term price diffusion depends on the limit order 
removal and this way is related to the market order flow. Furthermore, the empirical analysis 
confirms that the volatility calculated for the fixed number of transactions has the long 
memory properties as well and it is correlated with real time volatility ^| . We accumulate 
all these results into the assumption that standard deviation W(t, r d ) may be proportional 
to the square root of the trading activity, i.e., W(t,r d ) ~ ky/N(t, r d ). This enables us to 
propose a simple model of return 
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P cumulative 




FIG. 4: (a) Cumulative probability distribution function of the volatility, P > (v), averaged over 
10 intervals calculated from the series of n(t) generated by Eqs. (|12|) and ()18|) . all parameters 
are the same as in previous calculations. Dashed line approximates the power law P(v) ~ 1/v 2 ' 8 . 
(b) Power spectral density S(f) of v calculated from FFT of the same series n(t). Straight lines 
approximate power spectral density 5 ~ l//^ 1 - 2 with j3\ = 0.6 and /?2 = 0.24. 



x(t,r d ) ~kN(t,T A )e t (18) 

and related model of volatility v = \x(t, t^) \ based on the proposed model of trading activity 
(I12|) . We generate series of trade flow n(t) numerically solving Eq. (fT2|) with variable steps 
of time Atj = hi = n /ni and calculate the trading activity in subsequent time intervals 
as N(t, r d ) = j/ +Td n(t')dt' . This enables us to generate series of return x(t, r d ), of volatility 
v (t, Td) = \x(t, Td)| and of the averaged volatility v = — V {U, r d)- 
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In figure H] we demonstrate cumulative distribution of v and the power spectral density of 
v(t, Td) calculated from FFT. We see that proposed model enables us to catch up the main 
features of the volatility: the power law distribution with exponent 2.8 and power spectral 
density with two exponents (3\ = 0.6 and 02 = 0.24. This is in a good agreement with the 
empirical data jisl. fioj] . 

V. CONCLUSIONS 



Starting from the concept of the fractal point processes we proposed process driven by 
the nonlinear stochastic differential equation and based on the earlier introduced stochastic 



or events in the physical, biological and social systems when their statistics exhibit power- 
law scaling indicating that the represented phenomena contains clusters of events over all 
scales. First of all, we analyze the statistical properties of trading activity and waiting time 
in financial markets by the proposed Poisson process with the stochastic rate defined as 
a stand-alone stochastic variable. We consider the stochastic rate as continuous one and 
model it by the stochastic differential equation, exhibiting long-range memory properties 
[loj | . Further we propose a new form of the stochastic differential equation combining two 
powers of multiplicative noise: one responsible for the probability distribution function and 
another responsible for the power spectral density. The proposed new form of the continuous 
stochastic differential equation enabled us to reproduce the main statistical properties of the 
trading activity and waiting time, observable in the financial markets. In the new model the 
power spectral density with two different scaling exponents arise. This is in agreement with 
the empirical power spectrum of volatility and implies that the market behavior may be 
dependant on the level of activity. One can observe at least two stages in market behavior: 
calm and excited. Finally, we propose a very simple stochastic relation between trading 
activity and return to reproduce the statistical properties of volatility. This enabled us to 
model empirical distribution and long-range memory of volatility. 
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point process model 




may serve as a possible model of the flow of points 
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